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Abstract. I describe the current status of a collaboration with J. D. Romano, 
R. H. Price, and W. Krivan to model the geometry of and gravitational radiation 
emitted by a binary system of compact objects in the regime where non- 
perturbative gravitational effects exist, but the rate of inspiral is still small relative 
to the orbital frequency. The method of looking for a stationary spacetime 
which approximates the evolving solution is initially being tested on a simpler 
model with an additional translational symmetry. This report consists of a 
general description of the method, followed by summaries of three techniques in 
varying stages of development: the simplification of the Einstein equations in 
the presence of two commuting Killing vectors which form a non-orthogonally- 
transitive symmetry group, the boundary conditions appropriate to the balance of 
ingoing and outgoing radiation needed to reconcile a stationary radiating solution 
with conservation of energy, and the treatment of gravitational waves far from 
the sources as linearized perturbations to the Levi-Civita spacetime. The poster 



presentation with which this papei is associated is available on line at http://www 
itp.unibc.ch/~whclan /poster. ps. gz and the curre nt status of the project is described 



http: / /www-ftp. unibe.ch/ ~whclan/qsbi.html 



1. Introduction 

1.1. Inspiral of Compact Object Binaries 

A pair of compact objects (black holes or neutron stars) in binary orbit about one 
another is stable in Newtonian gravity. In general relativity, however, the system will 
emit gravitational radiation, causing the bodies to spiral in towards one another. The 
gravitational radiation given off by this system is a prime candidate for detection by 
upcoming gravitational wave telescopes such as VIRGO and LIGO. 

When the compact objects are far apart, all gravitational effects are weak and one 
can use the Post- Newtonian approximation, expanding in powers of GM/Rc 2 . The 
final "plunge", where the objects cease to orbit and collide rapidly, can be modelled 
by full general relativistic three-plus-one supercomputer simulations. An intermediate 
late inspiral phase, where strong gravitational effects are important, but the fraction 
of total energy lost to gravitational radiation each orbit is small, cannot be handled 
by supercomputer evolutions, which become unstable after several orbits. To handle 
strong gravity and slow inspiral, we need to use an approximation scheme. Modelling 
this intermediate phase will allow us to determine the early gravitational wave signal 
as well as provide initial data for the final plunge simulations. 
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1.2. Quasi- Stationary Approximation 

The idea, initially proposed by Steven Detweiler jj], is that if the inspiral is slow, the 
system is nearly periodic: after one orbit, the objects have returned almost to their 
original locations, and radiation which has moved out has been replaced with new 
radiation of approximately the same shape. If the objects' orbits are circular rather 
than elliptical, the spacetime is nearly stationary. If the approximate orbital frequency 
is f2, moving forward in time by St and rotating the resulting spatial slice by — f2 St 
will not change the picture very much. 

Our approach is to replace the true spacetime, which has this approximate 
continuous symmetry, with another solution to Einstein's equations in which the 
symmetry is exact. This "stationary quasi-solution" , which must somehow replace 
the energy lost in radiation to prevent inspiral, should approximate the physical 
"quasi-stationary solution" over some time interval. This replaces a three-plus-one- 
dimensional numerical evolution with a three-dimensional instantaneous solution, 
reducing greatly the computing power required and hopefully escaping some of the 
numerical instabilities associated with evolution. 

To further simplify the problem, we will initially look for a spacetime with an 
additional translational symmetry orthogonal to the orbital plane. This toy model of 
co-orbiting cosmic strings further reduces the numerical problem to two dimensions 
while hopefully retaining some of the qualitative features. 

The rest of this paper consists of summaries of three techniques being developed 
for this project, but which may prove useful in other settings as well. Section ^| 
describes research done in collaboration with Joseph D. Romano on the refinement 
of a formalism developed by Geroch || to simplify the Einstein equations when the 
spacetime admits two commuting Killing vectors, as is the case with the co-orbiting 
cosmic string model. Section || discusses the construction of a stationary radiating 
solution to a wave equation without the use of external forces by finding a preferred 
solution containing a balance of incoming and outgoing radiation, the subject of a 
collaboration with William Krivan and Richard H. Price Section [| briefly touches 
on the concept of describing small-amplitude gravitational radiation as a perturbation 
not to Minkowski spacetime but to the Levi-Civita solution || for a single cosmic 
string. 

2. QSBI I: Spacetimes with Two Killing Vectors 
(Whelan and Romano [^]) 

2.1. Lack of a Block- Diagonal Coordinate Basis 

The spacetime of two co-orbiting cosmic strings has two continuous symmetries, which 
are described by two Killing vector fields. One of these, which we call Kq, represents 
the combination of time translation and rotation about the orbital axis which leaves the 
spacetime unchanged. It can be thought of roughly in terms of traditional cylindrical 
coordinates as dt + Sidtp. The other Killing vector simply corresponds to translation 
along the strings and can be thought of as K\ ~ d z . 

In a numerical determination of the spacetime geometry, one seeks to fix the coor- 
dinate (i.e., gauge) information completely, and thus calculate the minimum number 
of quantities necessary to define the geometry. It is desirable, of course, to choose a 
gauge which takes advantage of the symmetries of the problem. Following from the 
example of static, axisymmetric spacetimes , we might wish to define two Killing co- 
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ordinates {x A } (x° ~ t, x 1 ~ z), supplemented by two other coordinates {x 1 } (x 2 ~ p, 
x 3 ~ tp ~ — fit) on a subspace orthogonal to the two Killing vectors, and bring the 
line element into a block-diagonal form 

X AB ({x k })dx A dx B + j ij ({x k })dx i dx j . (1) 

For our set of Killing symmetries, however, this fails because the symmetry group is 
not orthogonally transitive, which means that it is impossible to construct two-surfaces 
everywhere orthogonal to both Killing vector trajectories. The measure of this failure 
is a pair of scalar fields {ca} given by 

c A = *(K AK 1 AdK A ). (2) 

The co- rotational Killing vector Kq is not surface- forming (Kq A dKo ^ 0), which 
leads to the non- vanishing of cq, indicating a lack of orthogonal transitivity. 

2.2. Manifold of Killing Vector Orbits 

Although the metric components in any coordinate basis will have non- vanishing "cross 
terms" {gAi} between the Killing and non-Killing directions, it still possible to describe 
the geometry in terms of two matrices {Xab} and {"fij} along with the two scalars 
Co and c\. This done using a construction due to Geroch || which defines a two- 
manifold S of Killing vector orbits. The coordinates {x l \i — 2,3} are coordinates on 
this two-manifold, and any tensor on the original four-dimensional spacetime manifold 
A4 which has vanishing inner products with, and Lie derivatives along, the Killing 
vectors corresponds to a tensor on the two-manifold S. In particular, the tensors and 
scalars on S which describe the spacetime geometry are 

• The symmetric matrix of inner products Xab — Ka • Kb 

• The metric on S, which has components {7^ } and corresponds to the projection 
tensor g^ — X AB Ka^Kbv on M (where {X AB } is the inverse of {Xab}) 

• The two scalars ca which define the lack of orthogonal transitivity. 

Given these objects, it is possible to define a non-coordinate basis on the four- 
dimensional spacetime by eA = Ka and e\ = ^ijdx 3 . In this basis the metric has 
components 

gAB = Xab , gAi = , gij = 7^ ; (3) 
the {ca} give the commutation coefficients of the non-commuting basis vectors: 

[e 2 ,e 3 ] cx X AB c A K B . (4) 

2.3. Einstein Equations 

Geroch's original paper contained four partial differential equations involving {"fij}, 
Xoo, A01, An, Co, and ci which were equivalent to the vacuum Einstein equations 
G>„ = 0. We have found a streamlined derivation which treats the indexed quantities 
Xab and ca as single entities, rather than dealing with each component individually. 
We have also derived explicit expressions (given in ||]) for the components Gab, GAi, 
and Gij which can be used even when the stress-energy tensor is non- vanishing. 

In the case where the off-block-diagonal components {T^j} of the stress-energy 
tensor vanish, the Einstein equations GAi = simply say that the scalars ca are 



4 



constants. A convenient choice of gauge in which to solve the remaining six Einstein 
equations is a basis in which 

c EE 20 , Cl EE , (5) 

Aoo = (X + X(p, ip) 2 )Z(p, tp)- 1 , \ i=X(p,<p), \n=Z(p,<p), (6) 
722 = 1, 723 = 0, 733 = -X(p, Lp)~ 1 p 2 F(p,(p) . (7) 

These four functions X, Z, X, and F of the two coordinates x 2 = p and x 3 = ip 
provide a fixed gauge in which the problem can be solved numerically. There are six 
block-diagonal Einstein equations (for {Gab} and {G^}), but only four of them are 
independent because of the contracted Bianchi identities (dx l ) vG^ = 0. These 
equations will be combined with a treatment of the cosmic string sources and the 
boundary conditions at infinity and the origin to produce a numerically determined 
spacetime. 



3. QSBI II: Radiation-Balanced Boundary Conditions 
(Whelan, Krivan and Price [Q]) 

3.1. Radiative Boundary Conditions 

The true physical spacetime which we are ultimately modelling contains gravitational 
radiation at infinity which is outgoing. This loss of energy leads to decay of the orbits 
and inspiral of the compact objects in a non-stationary solution. In many radiative 
problems, a solution with outgoing radiation can be constrained to be stationary by an 
external force whose agent does not couple to the radiation. However, this is at odds 
with the idea that in General Relativity all matter gravitates, so we should look instead 
for a solution where there is no net energy loss to infinity due to the radiation. A naive 
replacement for the outgoing radiation boundary conditions would be a standing wave 
condition. However, that turns out to be inappropriate, as standing waves require 
a node (Dirichlet) or extremum (Neumann) at a particular location, and a standing 
wave condition thus fails to converge to a well-defined limit as that location is moved 
out to infinity. 

We have been investigating the question of how to implement a sensible boundary 
condition leading to a balance of radiation in the context of a simple theory: a 
nonlinear scalar field ip(t, p, (f>) in two-plus-one dimensions, where the source a and 
field ip are required to be co-rotating (i.e., depend only on p and ip = (j> — fit). This 
causes the field equation to take the form 

U^ = — +p - + (p -fic ) — = a + Xi,. (8) 

In the numerical solution for the non-linear equation (|8|), one uses a finite co- 
ordinate grid and specifies a set of boundary conditions at the large finite radius 



p = R which marks the end of the grid. Table 3.1 shows several possible boundary 
conditions. We call, e.g., the solution resulting from the application of the outgoing 
boundary condition at a particular R, i/' j U ^- The local conditions defining outgoing 
or ingoing radiation each produce well-defined limits (which we call simply ijj ou ^ and 

i/j lla ) as R is taken to infinity, but the standing wave solutions do not converge to any 
limit. 
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BCs at p = R 




R - 


— > oo limit 


Outgoing 


(dp + c^dt)^ 1 = 


^out 






Ingoing 


(d p - c-'dt)^ = o 






p—<p) ^ gimn(c _ V+*) 


Neumann SW 


d p i$ = 






N/A 


Dirichlet SW 








N/A 



Table 1. Some large-distance boundary conditions. Purely ingoing or outgoing 
radiation has a well-defined limit as the radius R at which it is applied is taken to 
infinity, but standing wave conditions (e.g., Neumann or Dirichlet) do not. 



Taking the average of the R — > oo limits ip oxa and defines a function 

^avg = i^in + ^outj (9) 

which has a balance of ingoing and outgoing radiation without reference to any 
particular radius R. In the linear (A = 0) theory, where the principle of superposition 
holds, it is a solution to the field equations, but in the non- linear theory it is not. It 
is also not defined as the limit of any local boundary conditions (in particular, the 
is average of out- and ingoing solutions is not produced by an average of out- and 
ingoing boundary conditions). 



3.2. Green's Function Solution 

The tasks of defining a solution analogous to ?/> av & and numerically determining that 
solution in the absence of a local boundary condition can both be accomplished by a 
Green's function method. 

If we define a Green's function G(p, ip\po, <fo) such that 

n 2 G = p Q 1 6(p-po)5(<p-<p Q ), (10) 

the differential equation (JsJ) can be converted into an integral equation 

^ip,v)= / PQdp d(p G(p,(p\pa,ipo) [a(p Q ,(p ) + \ip 3 {p,<p)] . (11) 



Solving to the Green's function equation (|10|) requires the application of boundary 
conditions at infinity, so there are actually a family of Green's functions with boundary 
conditions "built into" them. For example, using the retarded Green's function G ou * 
in ([ll]) gives the outgoing solution ^ out , using the advanced Green's function G m gives 
the ingoing solution ip 111 , and similarly for any of the boundary conditions defined at 



finite radii in Table B.l 



Since the Green's function equation (10) is linear, even when the wave equation 



(Q) is not, we can always average the retarded and advanced Green's functions to give 
a time-symmetric Green's function 

G sym = l (G in + G out ) (12) 

which defines a radiation-balanced solution ip s y m . 

In the linear theory, this solution ip s y m will be equal to the average ?/> av & of the 
ingoing and outgoing solutions. In the non-linear theory, we can calculate the two 
numerically (by iterating (fTy)) and compare them. We are ultimately trying to use 
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Figure 1. Comparison between i/) s y m and ^ av S in the wave zone for a source of 
two particles of opposite charge at p/2 = .25, if = {% rfc tt)/2, with a rotational 
frequency of Q = 1/2 and an orbital velocity of one-fourth the speed of light. The 
fields are plotted versus p/2 for the range of angles ip £ [l,7r/2]. The amplitude of 
i/) av S is around 20% smaller and the phase is shifted 40 degrees relative to ip s y m . 
This agreement is only middling because the sources are still somewhat relativistic. 



the outgoing piece of i/j^ 111 as an approximation for ip ou ^ (and thus as a further 
approximation for the actual inspiralling solution), and we expect that approximation 
to be good when -0 s y m w -0 av S. 



3.3. Numerical Results 

Numerical calculations for a strongly non-linear (A = 20) theory have shown that 
the approximation becomes good as the orbital velocity of the sources becomes non- 
relativistic. Figures [[jand^ show the comparison between ip s y m and ifj av & for orbital 
velocities of one-fourth and one-eighth the speed of light, respectively. 

This agreement occurs even when the theory is highly nonlinear because, even 
though ifj m and tp ou ^ differ in the wave zone, they are approximately the same in the 
inner, strong-field region, as illustrated in Figure ||. 

The final step in the scalar field theory analysis will be to numerically extract the 
"outgoing part" from the time-symmetric solution ■0 s y m and compare it directly to 
^out. we ex p ec t that the two will agree when ^> s y m is approximately equal to ?/; av S. 



4. QSBI III: Gravitational Waves on a Cosmic String Background 

4-.1. Schwarzschild vs. Levi-Civita Backgrounds 

In the three-plus-one binary inspiral problem, with non-extended sources, gravitational 
waves at infinity can be treated as perturbations to Minkowski spacetime. This is 
because at large distances, where higher multipoles of the non-radiative field can be 
ignored, the non-radiative part of the geometry is associated with the Schwarzschild 
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Figure 2. Comparison between ?p s y m and i/j av S in the wave zone for a source of two 
particles of opposite charge at p/2 = .125, ip = (7r±7r)/2, with a rotational frequency 
of Q = 1/2 and an orbital velocity of one-eighth the speed of light. The fields are 
plotted versus p/2 for the range of angles ip £ fl,7r/2]. With these less relativistic 
sources, the agreement is better than in Figure jj, with approximately a 9% smaller 
amplitude for i/j av S and 7 degree phase shift. 



metric 



(1 - r s /r) c 2 dt 2 + (1 - r s /r)- 1 dr 2 + r 2 dtf 



which contains a length scale r s 



2GM/c 2 ; taking r/r s 



(13) 

oo reduces the 

Schwarzschild line element to the Minkowski one, which is just another way of saying 
that the Schwarzschild solution is asymptotically flat. 

With infinitely extended sources, this is not possible. Far from the strings, where 
the internal structure is not felt, the non-radiative effects are approximated by the 
Levi-Civita solution Q of a single cosmic string at the origin with mass-per-unit 
length A, whose line element is given for small A by 



p 2C c 2 d? 



-ic 



dz 2 + dp 2 



(14) 



The mass parameter C — 2GA/c 2 is now dimensionless; this lack of a length scale 
means that Levi-Civita spacetime is not asymptotically flat and we must use Levi- 
Civita rather than Minkowski as our background metric. 



4-. 2. Perturbations to Levi-Civita 

Suppose the metric is given by the Levi-Civita metric plus a small perturbation 
hfj,v(t, p, 4>). (In our problem the t and (j> dependence will further be restricted to 
a function of — ilt, but we are thinking here of more general considerations.) We 
can divide the components {h^} into those which are odd under inversion of the z 
coordinate {h z t, h zp , and h zc j,) and those which are even (h zz , htt, h tp , ht^, h pp , h^, 
and h^c/,). In the linearized theory, the odd and even sectors are uncoupled. 
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Figure 3. Illustration of the agreement between jft ou ^ and i/> m in the strong- field 
region. The parameters are the same as in Figure bl i.e., oppositely-charged source 
particles at p/2 = .125, tp = (-7T ± i")/2, a rotational frequency of Q = 1/2 and an 
orbital velocity of one-eighth the speed of light. The fields are plotted against p/2 for 
ip = 3n/8. Note that while the two solutions differ in the outer, wave zone, they are 
nearly equal in the inner regions which are the source of the most of the nonlinear 
effects, leading to the approximate agreement of ip s y m and ?/> av § (which are identical 
in the linear theory) illustrated in Figure bL 



To count the physical degrees of freedom, we consider the effects of a gauge 
transformation 

h MV ->■ hft V - 2 (15) 

by a gauge parameter ^{t, p,<f>), as well as the Hamiltonian constraint H and the 
three components {H Zl H p , H^} of the momentum constraint. 

The three odd components of the metric perturbation consist of one gauge degree 
of freedom £ z , one constraint H z , and one physical degree of freedom, which can be 
thought of as the "cross polarization" (since it involves h Z( p). 

The seven even components include three gauge degrees of freedom £ t , £ p , and 
three constraints H, H p , and and one physical "plus polarization" involving h zz 
and h^. 

The next step in a treatment akin to that done by Moncrief for waves on 
a Schwarzschild background will be to derive explicit wave equations for these two 
gauge-invariant physical degrees of freedom. For the purposes of the quasi-stationary 
binary inspiral program, it may suffice to work with linearized versions of the gauge- 
fixed equations derived in [Q], but the picture of general perturbations to Levi-Civita 
will at least be useful in interpreting the quantities involved. 
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5. Future Outlook 

The work described in Section ^ has produced non-linear differential equations which 
can be numerically implemented given a suitable treatment of the sources and the 
boundary conditions on the radiation at the outer boundary. The study of radiation 
balance in the context of non-linear scalar field theory in Section || will be completed 
by the extraction of the outgoing part of the time-symmetric Green's function solution 
and the comparison between that and the actual outgoing solution (which we expect 
to be unavailable in the gravitational case). Application of these methods to General 
Relativity will be complicated by the fact that it is probably not possible to formulate 
a (linear) Green's Function solution to the non-linear wave equation in that case. 
The analysis of perturbations to Levi-Civita spacetime described in Section || is still 
in the early stages, and its full development may or may not be necessary to the 
understanding of our problem, which possesses the additional co-rotational Killing 
vector. 

Finally, should the technique of finding a stationary solution with a balance of 
ingoing and outgoing radiation prove successful in the essentially two-dimensional 
problem of orbiting cosmic strings, it can be applied to the computationally more 
intensive three-dimensional problem. 
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